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Abstract. Through a Euclidean path integral we establish that the density fluctuations of 
a Fermi fluid in one dimension are related to vicinal surfaces and to the stochastic dynamics 
of particles interacting through long range forces with inverse distance decay. In the surface 
picture one easily obtains the Haldane relation and identifies the scaling exponents governing 
the low energy, Luttinger liquid behavior. For the stochastic particle model we develop a 
hydrodynamic fluctuation theory, through which in some cases the large distance Gaussian 
fluctuations are proved nonperturbatively. 
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I. INTRODUCTION 



As pointed out by Haldane some time ago ]T|,0, spinless fermions in one dimension 



interacting through a short range potential have universal ground state correlations. The 



where low energy characteristics turn out to be labelled by two free parameters, traditionally 
denoted as the renormalized sound velocity v s and K. These parameters must be matched 
to the microscopic Fermi liquid under consideration. In fact, we will see that they are given 
by suitable second derivatives of the ground state energy/length. 

The two most prominent predictions of the Luttinger liquid scenario are 
(i) The momentum distribution behaves as 



close to the Fermi momentum hp- Compared to the noninteracting case, the Fermi fluid 
looses its gap at k F and the Fermi "surface" is retained as power law singularity only with 
anomalous exponent a = h(K + — 2). 

(ii) The density fluctuations in the ground state are severely suppressed and strongly corre- 
lated. Nevertheless they have Gaussian statistics. This is a consequence of the bosonization 
of the density field, which is the basic observation leading to the exact solution of the 
Tomanaga-Luttinger hamiltonian ||. 

The anomalous momentum distribution ( |1 . 1|) has been studied by G. Gallavotti and 
coworkers through a rigorous implementation of a renormalization group zooming onto 
the Fermi surface. The present paper focuses on the density fluctuations (ii). We recall 
first that for the ideal Fermi fluid the structure function, i.e. the Fourier transform of the 
density-density correlations is given by 



universal properties are computed on the basis of the Tomanaga-Luttinger hamiltonian 0-[5| 



(ar(k)a(k)) ~ \k — kp\ a sign(fc — k F ) + regular part 



(1.1) 




(1.2) 
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with kp = up the Fermi momentum and p the density. The interaction smoothens the cusp 



at 2kp, similarly to fll.lQ , as nas been proved recently by Benfatto and Mastropietro ||. 
For density fluctuations the behavior near k = is of interest. Here the interaction has a 
much less spectacular effect. It only changes the openening angle of the cusp (and of course 
modifies the straight piece). Here we plan to go beyond the static 2-point function and to 
study the small k, uo behavior of all n-point functions, including their frequency behavior, i.e. 
we plan to study the generating functional of the density field for large space-time distances, 
at which the Gaussian statistics should be recovered. 

As will be explained in detail below the density fluctuations are most conveniently in- 
vestigated through the path integral for the world lines of the fermions. This form leads to 
two other physical interpretations. One may think of the world lines as steps of a vicinal 
surface and use the statistical mechanics of surfaces. In this picture the Gaussian fluc- 
tuations are fairly immediate and Haldane's parameters are identified as suitable second 
derivatives of the surface tension. In the second interpretation one regards the fermionic 
world lines as trajectories of particles whose motion is then governed by certain stochastic 
differential equations. Such stochastic particle systems are usually described through a hy- 
drodynamic type fluctuation theory. In our case the forces between the particles decay like 
the inverse distance and are therefore long ranged. We will develop a suitable modification 
of the standard hydrodynamic fluctuation theory. In the framework of stochastic particle 
system, at least for some cases, we prove the Gaussian fluctuations without going through 
the perturbative double expansion in n and the interaction strength. 

Different physical interpretations of the same theoretical model lead to alternative ap- 
proximation schemes. Properties which look very deep in one formulation are physically 
obvious in another. We regard it as interesting that one- dimensional Fermi liquids allow for 
three distinct physical interpretations and try to explore their interconnections. 
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II. BASIC MODELS AND THEIR PATH INTEGRALS 



Let us start with the two prototypical models. 

(i) Fermions on a ring [0, £] interacting through a short range potential. The hamiltonian 
reads 

N ] 8 2 1 " 

^=-E^+2 £ V(xi-x s ). (2.1) 

3=1 1 i^j=l 

We have set m — 1 = K. V is a short range potential and Xj G [0, £] with periodic boundary 
conditions. By the Pauli exclusion principle, the ground state wave function, has to vanish 
at {xi = Xj} (Dirichlet boundary conditions). For density fluctuations the sign changes in 
ijj do not matter. Therefore, equivalently, we may regard Q2.1|) as the hamiltonian of bosons 
with hard core. Formally, this corresponds to adding an infinitely strong repulsive 5-potential 
to V. 

(ii) Fermions on the periodic lattice [1, . . .,£]. In second quantized form the hamiltonian 

reads 

i 

H = ^{-ala x+1 - al +1 a x - Aala x al +1 a x+1 } , (2.2) 

x=l 

a i+i — a i- For future use, we stated only the particular case of a nearest neighbor interaction. 
Through the Jordan- Wigner transformation (|2.2|) turns into the XXZ hamiltonian as 

1 1 

H = ^{-2^*1+1 + <&Ui) ~ (2-3) 

x=l 

with periodic boundary conditions at+i = where a x are the Pauli spin matrices at site 
x. 

Both hamiltonians, ( |2.1| ) and ( |2.3| ), generate a path integral for the the statistical weight 
of fermionic world lines. For ( |2.1| ) the free measure, P Q , are N independent Brownian motions 
Xi(t), . . . , x n (t). By the Pauli exclusion principle they are constrained not to cross, i.e. 

Xj(t) < x j+1 (t) mod i for all* , j = 1, . . . , N, (2.4) 
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which ensures the Dirichlet boundary condition at {x{ = Xj}. The statistical weight of the 
world lines is then given by 

1 1 N f 

-P ex V [--J2 dtV( Xi (t) - Xj (t))]xNC , (2.5) 
¥i=i 

where Z is the normalizing partition function and Xnc is the indicator function of the set 
defined in (|2.4f) . xnc restricts the path integral Q2.5|) to non-crossing paths only. In the 



bosonic language this constraint corresponds to an infinitely repulsive 5-interaction. 

The statistical weight generated by ( |2.3|) is the same as in (|2.5|) . Only the Brownian 
motion Xj(t) is to be replaced by a continuous time random walk on [1, ... , l\ with jump rate 
1 to the right and to the left neighbor. In our example we picked the particular interaction 
potential V(x) = —A for \x\ = 1 and V(x) = otherwise. 

One physical picture suggests itself immediately. We can think of the fermionic world 
lines as the trajectories of a stochastic particle system. In fact, in the limit t — > oo it will 
be a stationary Markov process. It is a diffusion process in case (i) and a jump process 
in case (ii). The quantum mechanical free Bose field of density fluctuations corresponds 
to space-time Gaussian fluctuations in the stochastic particle system. This is the usual 
hydrodynamic fluctuation theory as governed by a linear Langevin equation. Unfortunately, 
the non-crossing constraint results in repulsive long range (l/ic)-type forces between the 
particles. Therefore the standard hydrodynamic theory, developed for short range forces, 
does not apply. We will explain in Section [TV] the required modifications which are in fact 
surprisingly little. 

An alternative physical picture is also well known PJITfl, but less immediate. We think 
of Xj(t) as a step of unit height for a vicinal surface. More explicitely we introduce a height 
function h(x, t) for the height of a crystal surface relative to the x, t reference plane. Then 

8 N 

—h(x,t)=Y,S(x-x 1 (t)). (2.6) 

i=i 

In the t-direction the average slope vanishes, whereas in the x-direction it is given by the 
particle density p = N/£ . Periodic boundary conditions for the particles means that the 



surface is extended by repetition at average slope N/i along the 

The stepped surface h(x, t) fluctuates with a statistical weight given by ( |2.5| ). In general, 
surface fluctuations are expected to be governed by a free massless Gaussian field in the 
infrared limit with a strength determined by the matrix of second derivatives of the surface 
free energy. Since it coincides with the ground state energy of the hamiltonian ( |2.1| ), we 
have a direct way to identify the parameters in the Luttinger hamiltonian. The Haldane 
relation |l| then follows simple consequence. 

The surface picture also indicates the limitations of the Luttinger liquid concept. The 
surface free energy may have cusps and/or flat pieces. The former case corresponds to a 
roughening transition where logarithmic fluctuations are suppressed to order one fluctua- 
tions. The latter case is step bunching. The slope "segregates" into macroscopic regions. 
The steps are closer together than expected by naive counting. 

To give a brief outline. In Section |T| we develop the surface picture in more detail. In 
Section [TV] we explore the hydrodynamic fluctuation theory for the dynamics of world lines 
and its consistency with the surface picture. We argue for the validity of a linear Langerin 
equation governing the density fluctuations. The approximations can be controlled for the 
Calogero-Sutherland model and for a general system with short range interaction provided 
the two-point function scales. The necessary computations are provided in an Appendix. 



III. SURFACE FLUCTUATIONS AND BOSONIZATION 

For the sake of concreteness we first discuss fermions on a lattice, cf. hamiltonians ( |2.2|) 
and ( |2.3| ). The relations derived below are general, however. We use the <r 3 -representation 



and it is convenient to set up a corresponding notation. We define rj x = (1 + cr|)/2. Then 
7] x = 0, 1 and we interpret r) x = 1 as a surface step at site x. A whole step configuration is 
denoted by rj. Then H of ( |2.3| ) in the ^-representation becomes a linear operator acting 
on functions f(rj) and is given by 



5 



Hf(v) = ~ J>* - V* + i) 2 f(v xx+1 ) ~ A »fe Wfo) , (3-1) 

x=l x=l 

where the periodic boundary condition 7^ +1 = t]i is understood. rj xx+1 stands for the config- 
uration 7] with occupations at sites x and x + 1 interchanged. The path integral is generated 
by the transfer matrix {e~ tH ) vv ' , t > 0. For an isolated step the first term in ( |3.1| ) yields 
a symmetric, time continuous random walk with jump rate 1 to its neighbors. For several 
steps (i] x — r/ x+ i) 2 ensures the non-crossing constraint. —A Yl x =i VxVx+i is a potential. Thus 
we see that ( |3.1| ) indeed generates a path integral of the form ( |2.5| ) with Xj(t) replaced by 
rj x (t) and Pq standing for independent random walks on the lattice [1, . . .£]. For later use 
we define the steps in space-time by 

rj x (t) = 1 (0), if at "time" t there is a (no) step at site x, 

x = 1, . . . ,£ and < t < T. The steps have a statistical weight given by ( |2.5|) . Clearly 
the number of steps is N = Y^ x =i Vx(t) independent of t. In the dynamic picture we regard 
particle and r\ x = as no particle at site x . The surface steps are then the 
world lines of the particles. Particles jump to their neighboring sites, but are never created 
and destroyed. We will use 'step' and 'particle' interchangeably. 

In the crystallographic literature our model is known as terrace-step-kink (TSK) model. 
It describes a high symmetry crystal surface miscut by a small angle. Such a vicinal surface 
consists of a regular array of monoatomic steps. Through thermal activation the steps 
meander but they do not cross or terminate. The slope of the vicinal surface imposes the 
step density and their average orientation. The terraces are the constant height pieces 
between steps and kink refers to a step corner. 

The surface defined through ( fOl) , with Xj(t) replaced by f] x {t), is tilted in the ^-direction 
with slope p = Nj I. A complete picture emerges only if we tilt the surface also along the t- 
direction. To do so, we define the step current J xx+ i(t) through the bond (x,x+ 1). J xx+ i(t) 
is a sequence of (^-functions located at those times when a step jumps between x and x + 1. 
The 5-function carries a weight +1 (—1) if the jump is from x to x + 1 (x + 1 to x). The 
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tilt along the t-coordinate is enforced by the additional constraint 

t T 

£ [ Jxx+i(t)dt = NaT , (3.2) 

which implies that on the average each step has the drift a. 

If the step variables r] x (t), x — 1, . . . ,£, < t < T, are given, then by definition 

h(x + l,t) -h(x,t) =%(*). (3.3) 

731) can be integrated to yield 



X 

h(x,t) = ^^(0) - / J xx+ i{s)ds 

v=l I 



(3.4) 



y= L o 

with an arbitrary choice for the constant of integration. 



At this point it is convenient to go from the canonical to the grand canonical prescrip- 
tion. The variable conjugate to the number of particles is the chemical potential fi and we 
introduce A as conjugate variable to the total current (|3.2| ). Thereby are obtain the transfer 
matrix {e~ tH ) m ' with hamiltonian 



Hf(v) = ~ E^^ 1 - + e ~ A ( 1 - Vx)Vx + i)f(v xx+1 ) + &VxVx + if(v) ~ Wxfiv)} 

x=l 

— H t {n, \)f{r)) . (3.5) 

The exponential of the current gives weight e A to a right jump and weight e~ x to a left 
jump and the exponential of the particle number yields the potential fiJ2rj x . ( ^.5| ) is the 

X 

hamiltonian of the asymmetric XXZ model. 

To make H symmetric we have to continue A analytically to the imaginary axis, i.e. to 

replace A by — iX with A real. Tracing back to the Fermi hamiltonian ( |2.2|) we obtain 

e 

H = ^{-e'^a^i - e~ lX al +1 a x - Aa\a x a} x+l a x+ i + ^a\a x } . (3.6) 

x=l 

This means that the dispersion relation — 2cosfc is replaced by — 2 cos(A; — A). In the low 
energy limit the dominant contributions have a a total momentum (N)X. Thus the (ana- 
lytically continued) A regulates the average fermionic current on the ring [1, . . . ,£], X = 
corresponding to zero current. 



The thermodynamics of the surface is governed by the surface tension er(u) depending on 
the slope u = (p, pa), a is convex up. To relate it to the hamiltonian ( |3.5| ) it is convenient 
to define the Legendre transform, a, of a by 

<r(v) = inf(o-(u) - v • u) . (3.7) 

u 

a is convex down and in terms of H in ( |3.5Q it is defined by 

Hf*, A) = - Jim trexp[-T^(p, A)] . (3.8) 

If we take first the limit T — > oo, then ( |3.8|) becomes Ee(/j,, X)/£ with the ground state 
energy of if^ . Thus the surface tension is simply the Legendre transform of the ground 
state energy of Hg per site, 

<7(p,A)= lim ]Et(ji,X). (3.9) 

Thermodynamic fluctuation theory suggests that a small height fluctuation, 5h, relative 

to the average slope u has a probability proportional to 

I 2 

exp[-- ^(n)(ViSh)(V 3 5h)] , (3.10) 

where ^-(u) = <9V(u) / 'duiduj. Thus on a large scale the height fluctuations should be 
Gaussian with a covariance 

2 

(Xl^iMW 1 (3.11) 



in Fourier space. (|3.11 ) is the covariance of a free massless bosonic field in its Euclidean 



version. Fluctuations of two-dimensional surfaces grow logarithmically and are therefore not 
stationary. They become stationary by taking an x-derivate, which according to (|3.3j ) yields 
the step density. In Fourier space we only have to multiply by k\. The thermodynamic 
fluctuation theory for surfaces predicts the low energy behavior of the density fluctuations 
for the world lines of ( |3.5| ). 

We find it convenient to keep the t-dependence and use k for the Fourier transform with 
respect to x. Then, according to ( [3.1 1|) and (|3.3j ), we have 
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(Vx{t)Vx'{0)) ~P 2 ^ I dke^-^c^e-^- 1 ^ 1 (3.12) 



2 

for large \x — x'\ and \t\, where the parameters are defined by 

en = (y 2 + v 2 )/yc , 0"i2 = -v/yc , <7 22 = I/7C (3.13) 



The density fluctuations are Gaussian on the scale where fl3.12|) is valid. 



For the expression in (|3~T0|) to make sense cr(u) must be twice differentiable at u and the 
matrix of second derivatives D 2 a(u) > 0. Already for the simple model (3.5) with nearest 



neighbor step-step interaction only, this condition is not always satisfied, a is known from 
the Bethe ansatz [11|. For A > —2 at p = |,a = 0, the steps align in antiferromagnetic 



order. Thereby surface fluctuations are strongly suppressed, from logarithmic order to order 
one. Changing either p or a destroys this roughening transition. We refer to [|nj for the 
behavior close to the transition. 

On the attractive side, A > 0, steps may bunch to give macroscopic patches of slope 
u = (1, 0) and of slope u = (0, 0). This phase is bordered by the stochastic line, where —H 
of ( |3.5| ) generates a stochastic time evolution |TT[. Clearly, the condition is 

A = 2coshA (3.14) 



for all p. To reexpress Q3.14 ) in terms of a, we need a = de/dX, which is not known in 



closed form, however. For small A from linear response we have a = p(l — p)2sinhA, which 
implies that for small a 

A c = 2 + (a/2p(l-p)) 2 . (3.15) 



Coming back to ( |3.12|) we note that the density fluctuations decay as y\k\ and propagate 



with velocity v. The static (t = 0) covariance is \c\k\ which reflects that steps are stiffly 
arranged because of the non-crossing constraint. The parameters y,v,c can be expressed 
through the ground state energy as 

c 2 = det/^ = (d 2 e/dp 2 )(d 2 e/d\ 2 ) - (d 2 e/dpd\) 2 , 

7 = c /(d 2 e/dp 2 ) , v = (d 2 e/dpdX)/(d 2 e/dp 2 ) , (3.16) 



where we used (|3.9|) , ( |3.13| ). These relations are valid in general. 

Haldane [1|] observed that for Luttinger fluids the parameters in the low energy effective 



bosonic action are not independent. This is easily rederived from ( 3.16 ). In the notation of 
Haldane 

v s = 7 , v N = 5/i/Sp . (3.17) 
The average current is j = de/dX and in linear response j(X) = j(Xo) + Vj(Xq)(X — Ao), i.e. 

Vj = d 2 e/dX 2 . (3.18) 

The Haldane relation reads 

v N vj = v 2 s , (3.19) 

i.e. 

{d 2 e/d^)-\d 2 e/dX 2 ) = 7 2 , (3.20) 

which is to be compared with 

(d 2 e/dfi 2 )' 1 (d 2 e/dX 2 ) =1 2 + v 2 (3.21) 

by (|3.16|) . Thus the Haldane relation holds at v — 0, equivalently at A = 0, which he used 
implicitely by setting j(Xo) = in the linear response. The parameter K mentioned in the 
introduction is given by 

d 2 e d 2 e^_ 1/ 



K = = 7^)" V2 ( 3 - 22 ) 



For A = 0, the Hamiltonian ( |3.5| ) is symmetric. In the dynamic picture we have a time- 
reversible jump process for the particles for which detailed balance is satisfied. As in the 
short range case this gives rise to an Einstein relation and the Haldane relation can be 
viewed particular case. 
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For Brownian steps, as governed by the Hamiltonian ( [2.1| ) with Dirichlet conditions at 
{xi = Xj}, some simplifications compared to the general case occur. The drift (= tilt along 
t) is enforced by the constraint 



dt Xj (t) = aT, (3.23) 



o 



j = 1, . . . , N. Such a drift can be trivially removed by the global change of coordinates 
Vj(t) = Xj(t) — at. In contrast to the terrace-step-kink model, the tilting costs only in 
elastic energy for each step individually. Let E N be the ground state energy of ( |2.1| ) and 
define 

e(p)= lim -E N . (3.24) 

N,e^oo,N/e= P £ 

Then the ground state energy per unit volume at drift a is given by 

e(p, a) = e(p) + ^pa 2 . (3.25) 



Comparing with ( |3.13|) and using u = (p, ap) we obtain v = a, as expected, and 



c = ^p/e"{p) , 7 = vV(p) • (3-26) 

The Haldane relation is cy = p . 

We conclude this section by stating a precise conjecture. For the terrace-step-kink model 
we average over e~ l sites with some smooth test function /. From the surface picture, the 
x and t coordinates must be on equal footing. Hence, time is also scaled as e~ l and we 
introduce the fluctuation field 

C(f,t) = J2f(^)(^ 1 t)-p). (3.27) 

X 

Note that we are not in the standard central limit theorem setting. We sum over e~ l sites, 
but expect a fluctuation of order 1 only. In the same spirit, for the Brownian steps we define 
the fluctuation field as 
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£ £ (/^) =Y l f{ex j (e-H))-p J dxf(ex) . (3.28) 



Conjecture. We consider the path measure generated by the hamiltonian ( \3. 1\ ) in the limit 
T — > oo at fixed tilt a, cf. ( [<?. 2\ ), and in the limit I — > oo at fixed density p = N/£. This path 



measure governs the fluctuation field ( \3.21\) . Ife(p,a), resp. e(p,X), is sufficiently smooth 
(at least twice differentiable) , then 

lim r(/,*)=e(/,*) (3-29) 
in distribution. The limit is jointly Gaussian with covariance 

(£(/,*)£M)> = f dk l -c\k\ e-^' 4 l e- M f(k)*g(k). (3.30) 



The parameters c,^,v are given by (\3.16 ). Correspondingly, for fermions in the continuum 



with hamiltonian ( \2.1\ ) and path measure ( \2.f\ ), the fluctuation field ( 3.2fy) satisfies the limit 



fljOg j with v = a in (fiHO) . 



IV. HYDRODYNAMIC FLUCTUATION THEORY 

In the limit T — ► oo, for fixed £, the statistics of the steps become stationary in t and 
are governed by a stochastic process which by construction is Markov. To determine its 
generator let ip be the ground state and E the ground state energy for H of either (|2.1| ) or 
( |3.1|) , which satisfy 



Hip = Eip. (4.1) 

In case of ( |2.1|) we have if) > except at coinciding points Xj = x iy i ^ j, where i/j(x) = 0. 
For ( |3.1| ) we first have to specify the sector Y^ x =i = N. Then, for fixed N, e~ tH has a 
strictly positive integral kernel and by the Perron-Frobenius theorem the ground state ip is 
unique and ip > 0. The (backwards) generator of the Markov process for the steps is defined 
by 

Lf = -r\H-E)^f (4.2) 
12 



as acting on functions / over the configuration space. The Markov process has ip 2 as unique 
invariant measure. 

We carry out this construction for the XXZ model in the a 3 representation and obtain, 

in the notation of fl3.1| ), 
i 

Lf(v) = - + e ~^ 1 " % if(v xx+1 ) ~ m\ ■ ( 4 -3) 

x=l 

L is the generator for a stochastic lattice gas, where particles jump to their neighboring 
sites. The exchange rate, c xx+ i(rj), between sites x and x + 1 is given by 

c xx+l ( V ) = ^ xx+1 )/m- (4-4) 

The jumps to the right are biased by the factor e A and those to the left by e _A . If A = 0, H 
is symmetric and the stochastic evolution satisfies detailed balance. 

The rates ( |4.4j ) are determined through the ground state, which is not known in general. 
Only for A = one has the explicit ground state if). If we denote by x%, . . . , x^ the positions 
of the particles in the sector Ylt=i Vx = N, then in this sector 

N 

if>(x!, . . . ,x N ) = ] [ \mi(n(xi-Xj)/i)\. (4.5) 
i<j=i 

Therefore the exchange rates are given by 

(N) r / x \ — ^ .sm(n(y — x + 1) /i) . 

c x N x {i = eM'(Vx + ^Vx) E (bg( y^ {3/ _ x) fi ) ' ))%]■ (4-6) 

Taking formally the infinite volume limit yields 

c xx+1 = exp[-(r) x+ i - r) x ) Y] {log{- — ^—))r) y \ . (4.7) 

A — ' y — x 

y^x,x+l 

(|4.7|) teaches us several points. The rates are repulsive: more particles to the left of x 
favors a right jump of the particle at x. The rates are long ranged, which means that the 
finite range intuition is no longer applicable. Since log((y — x — l)/y — x) = —{y — x)~ x , at 
infinite volume the rates may be infinite. The dynamics with rates ( j4.7| ) is then defined only 
for |^| 2 almost all configurations. Because the case A = maps to a free fermion theory, one 
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can construct the infinite volume ground state as a measure on particle configurations and 
the Markov semigroup e Lt in the corresponding L 2 -space [|I5L[Ij|. This implies the almost 



sure existence of the dynamics. If one adds the nearest neighbor interaction A, one expects 
exchange rates to be qualitatively similar to ( |4.7| ). We are not aware of any result in this 
direction. 

For the Brownian steps L is the generator of a diffusion process given by 

N f) f) 2 

L/(*)=5>(*)— + (4.8) 

j=i i j 

x = (xi, . . . , Xn). The drift on the j-th. particle is 

ai = -^-logV. (4.9) 

Thus log ip is the potential for the diffusion process. 

For Brownian steps with exclusion only, i.e. V = 0, one finds that 

N 

CLj(x) = ^ cot(vr(xj - Xi)/£) . (4.10) 
For I — > oo this expression converges to 



i.e. to a repulsive 1/x-force as for the lattice gas. The model with drift (|4.11|) was introduced 
by Dyson Its potential is formally given by 



-^2log\xi-Xj\. (4.12) 



In fact, Dyson adds a confining external potential, which is quadratic and considers ( |4.11| ) 
for finite N. 

We conclude that density fluctuations in the fermionic ground state may equally well be 
studied through the stochastic dynamics (|4.3|) , (|4.8|) . If the potential is short range, there 
is a well developed machinery known as hydrodynamic fluctuation theory. One argues (and 
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proves in many model systems fl6|J17|| ) that density fluctuations relax diffusively and are 
driven by a white noise random flux. This corresponds to the linear Langevin equation 

= a J^cm) + ^ w ^t) . (4.i3) 

D is the bulk diffusion coefficient and o is the bulk mobility. They are related through 
a = \D with x the static compressibility of the lattice gas. W(x,t) is space-time white 
noise. The spatial derivative ensures the conservation of the number of particles. The 
stationary measure for ( |4.13|) is Gaussian with covariance matrix {a/ 2D). In passing we 



mention that in the short range case fl4.13|) holds only for A = 0, i.e. in the symmetric case. 
For A 7^ one switches over to the KPZ universality class. It is characterized by the dynamic 
exponent z = 3/2, rather than z = 2 as in ( |4.13|) . The fluctuations are non-Gaussian [|i~9f . 

Returning to the long range case of interest here, we have already obtained the limiting 
Gaussian fluctuations in the Conjecture. We write £(/, t) = J dxf(x)^(x,t) and note that 
( p.30| ) determines a semigroup in t. Therefore ( |3.30| ) must be the (stationary) solution of 



the linear Langevin equation 



J^(x, t) = (- 7v /-97&r 2 - v-j^)t(x, t) + Vn-^W(x, t) , (4.14) 



which is a surprisingly minimal modification compared to the short range case (|4.13 ). 7 
plays the role of D and c the one of a. The crucial difference is that a Fourier mode e tkx 
decays as e~ 7 ' fc "*' rather than as e~ Dk2 ^. 



V. SCALING LIMIT 



We support our general conjecture by arguing that the stochastic particle evolution is 



close to the Langevin equation ( |4.14| ). We will do so on a fairly formal level. In particular, 
we simply work in infinite volume. In the appendix we explain how parts of our arguments 
can be made rigorous. 

Our strategy is most easily explained for the Brownian step model. The equations of 
motion are 
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^x j (t) = a j (x(t)) + b j (t), (5.1) 

where bj(t) is normalized white noise independent for each j. The scaled fluctuation field 
satisfies then the differential equations 

j t C(M = J2f(x%t)) aj (x%t)) +eJ2f"(x%t)) + yfeY^f\x £ m h i(t) ■ (5-2) 

3 3 i 

Here x e At) = exj(e~H) and we used the scale invariance of white noise as bj(e~ l t) = \/sbj(t). 

The second term in (|5.2| ) converges to p j dxf"(x) = 0. The third term converges to a 
space-time Gaussian measure with covariance 

5{t-t')p [ dxf'(x)g'(x). (5.3) 



This uses only that with respect to the distribution given by \ip\ 2 we have f(exj) — > 

3 

p j dxf(x) in probability as e — > 0. (|5.3|) is in accordance with ( 4.14j) , since cy = p. Thus 

we are "only" left with to show that (recall that v — 0) 
t t 
J dsJ2f'(x s j {s))a j {x £ {s))ds~-J J2^V-d 2 /dx 2 f(x £ j (s))ds. (5.4) 
o o J 

If so, the integrated version of ( |5.2[ ) becomes a closed equation and agrees with ( |4. 14| ) . 



To establish ( p.4[ ) is certainly the hard part of the matter. In Appendix A we prove that, 
if we assume the 2-point function to scale, the substitution ( |5.4|) holds. This is of interest 
because the scaling of the 2-point function by itself implies already that the fluctuations are 
Gaussian. I would not know how to obtain such a result otherwise. 

A further example for wich the substitution ( |5.4|) holds is the Calogero-Sutherland model 
18| , where the pair potential is proportional to 1/x 2 . If the particles are on a ring [0, £] and 



if we take all the image potentials into account, then the ground state wave function of the 
Calogero-Sutherland model is 

N 

*{j(x 1 ,...,x n )= 11 Ismiirix.-x,)/^ 2 , (5.5) 
i<j=i 



(3 > 0. By (|4.9D the corresponding drift is given by 
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(3 N 71 

a j {x u ...,x N ) = - ^2 j cot(ir fa - Xj)/E) (5.6) 



2 



which in the limit I — > oo yields 



^ ^ ^^^^ O / > ■ . / > ■ 

Compared to (|4.8| ) only the strength is changed. 



^01 



The ground state energy per length is e(p) = hf3 ix p (note that our kinetic energy is 



-±A). Therefore, by ( ggp 



c=— l = \*Pp. (5.8) 
7rp 2 



For the Calogero-Sutherland model the substitution in ( |5.4| ) holds without averaging in 
time. This is somewhat surprising and very particular for the 1/x 2 potential. To complete 



the argument one needs the central limit result of Johansson |20f . 

Urn <exp[£ / (<■*,-)]> = exp[ii J dk\k\\f k \ 2 } (5.9) 

for smooth test functions with j dxf(x) = 0, where (•) is the average in the ground state 
( |5.5p in the infinite volume limit. The instructive computation is explained in Appendix B. 
Turning to the terrace-step-kink model the situation is more complicated, as may be 



anticipated from the short range exchange rates pl| . We consider the symmetric case and 
set A = 0. Using the method of martingales one can show that the variance of the time- 
integrated noise is given by 

Mv))t = (^(rMr/ 01 )^-^) 2 )*, (5-10) 
v 

which is to be compared with the prediction 

° 1 = dX 2 ~ |a=0 (5 ' n) 



from ( f4.14| ). Note that v = 0, since A = 0. By second order perturbation theory in A we 
obtain 
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d\ 2 

V 



- E E ^ 01 )(Vo - mXH - E)-\ Vx - Vx+1 )H V M ) . (5.12) 



x r\ 



Thus (|5.11|) holds only if the second term vanishes. 



For A = the sum YliVx — Vx+i)ip{v xx+1 ) is the total current J acting on if). Since 

x 

J if) = 0, we conclude that for free fermions (|5.11| ) holds. In this case the two-point is 
explicitely known and we can use the argument of Appendix A to prove the Conjecture. 
Alternatively the n-point density correlations can be written in terms of 2-point functions. 
By applying the closed loop theorem one again concludes that the scaling limit is Gaussian 
with covariance ( [PUD 

If A 0, the second term in (|5.12 ) is not expected to vanish. This can be verified by 



expanding (H — E)^ 1 to second order in A. Thus the martingale term ( |5. 1U| ) yields the 
wrong noise strength in the Langevin equation - the hallmark of the so-called non-gradient 
systems. The drift term 

t 

JdsJ2 f Mc^+iM^ 1 *)) (5-13) 
o x 

can no longer substituted deterministically as in ( |5.4j ). Somewhere hidden there must be an 
extra noise term. While for short range lattice gases this mechanism is understood by the 



beautiful work of Varadhan and Yau [21], cf. also [23], the situation looks rather complicated 



for the long range case considered here. 



VI. CONCLUSIONS 

We tied together three, at first sight, disconnected pieces: the Luttinger liquid behavior at 
low energy, surface fluctuations, and the hydrodynamic fluctuation theory for the stochastic 
dynamics of world lines. 

In the surface picture one can easily identify the universal low energy limit of the density 
fluctuations for one-dimensional Fermi fluids. Of course, as an input one needs a variant of 
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the Einstein fluctuation formula. The remainder of the argument is then straightforward. 
In particular, we show that the parameters of the Tomanaga-Luttinger Hamiltonian must 



be matched to suitable second derivatives of the energy/length, c.f. Eqs. ( |3.16| ). To our 
knowledge, this has not been discussed with sufficient clarity before. 

In the interpretation of the world lines as stochastic dynamics the particles have long 
range interactions, which result from the Dirichlet boundary condition at coinciding positions 
and thus from the Pauli exclusion principle. We develop a novel fluctuation theory for 
the long range case. In fact, the only modification is to replace in the drift term the 
Laplacian, —d 2 /dx 2 , by the nonlocal integral operator \J— d 2 /dx 2 . Heuristically, following 
the arguments of Tomanaga, the \k\ decay comes from linearizing the dispersion relation of 
the particles at the Fermi surface. The not so obvious point is that an interaction changes 
only the prefactor 7 but not the decay law itself. At least for the Calogero-Sutherland 
hamiltonian, the mechanism behind such a renormalization could be fully elucidated. 

APPENDIX A: 

For the Brownian steps we define the structure function, S(k,t), by 

<(£/(*i(*)))(Z>(*i))> = / dkf(ky g (k)s(k,t) (A.i) 

with j dxf(x) = = j dxg(x). The average is in the stationary process at infinite volume 
with density p. The existence of this limit and its spatial ergodicity is assumed here. In 
principle, this could be avoided by considering a finite circle of length £, for which the 
stationary process exists, and by scaling £ = e -1 /o -1 ,iV = e~ x . We adopt this strategy 
for the Calogero-Sutherland model in Appendix B. In essence, Fourier space becomes then 
discrete, k 6 27rZ. 

Our real assumption is the scaling of the structure function as 

YmLE- x S(ek,e-H) = -\k\e~^ km (A.2) 

e^O 2 



with c = \J p/e", 7 = y/pe". We want to prove that (|A.2|) implies 
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lim([ J dsj^aj^is)) f(xp)) + J ds £ 1 ^¥Jdtf f{x*{s))f) = 0. (A.3) 

•? J 

We note that the time average is needed. This reflects that the system takes a while to 
adjust to local perturbations. 

We work out the square and use that L is symmetric, (F(LG)) = (G(LF)), and 



"(/)(*) = £/(*;)■ ( A -4) 



Under our scaling the term n(f") vanishes. We have 

t t 

d Sl J ds 2 (f(f, Sl )f(f,s 2 )) 



t si t 

= J d Sl ( J ds^fine^-^fif)) + J ds 2 (f{f)e L ^- s ^f{f))) 

si 

t si t 

= J d Sl (- J ds 2 ^L(f(f')e L ^- s ^- 1 n%f)) + J ds^ifiDe^-^- 1 ^ 

si 

t 

= -2t(f(f'X(f))+2 J ds^- s (rf{f)e Ls ^rf{f)) 

o 

= -2t(f(f'X(f)} + 2(n £ (/, tX(f)} - 2{rf{f)rf{f)) . (A.5) 
For the first summand we have 

d _ d 



*«(/)» 



(f(f>V)) = ~e- l \ E(^ £ (/)^(/)) • ( A - 6 ) 

j j j 

Therefore in the limit e^Owe obtain 

J dk\f\ 2 (t\k\ 2 p - c\k\(l - e-^H*l)) . (A.7) 

The second term reads, with g = —^^—d 2 /dx 2 f 1 
t t 

-2 J d Sl J ds 2 (n%g, Sl )f(f',s 2 )) 



o o 
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t si 



-2 j ds 1 [Jds 2 (n £ (g)e L ^- s ^ 1 f(f)}+ J ds 2 {f(f)e L ^- s ^ n%)>] 



si 

t si t 



= -2 J ds x [- J ds 2 ^-(n%g)e L ^ s ^- 1 n%f)) + J ds 2 ^-(f(f')e L ^- s ^ n e (g))] 



si 



= -2 y ^[-2<n%K(/)> + 2<n e G/)e L - e ~V(/)>] . (A.8) 



o 

By assumption this expression converges to 



J dk\f\ 2 (-2t 1C k 2 + 2\k\c(l - e-^)) . (A.9) 

The third term is given by 

t t 

J ds! J ds 2 (n £ (g,s 1 )n £ (g,s 2 )) , (A.10) 





which by assumption converges to 

t t 



dSl / dS2 /*|/|V>|V* 1 ..-^/*|/r (( ^-c|,|(l-e--)). (A,!) 





Using that 7c = p, the sum of the three terms vanishes, as claimed. 
To complete the argument one notes that 

t 



3 

is a martingale with square 



M £ (f,t) = J dsJ2f(^-^))dbj(s) (A.12) 



t 

M(f,t) 2 = f dseJ2f(^A^)) 2 + M!(f,t), (A.13) 
i 

where M[(f, t) is again a martingale. We assume now the validity of a law of large numbers 

as 

l™e^0(exj) = P J dxg(x). (A.U) 
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Then, under suitable tightness, in the limit e — > 

M(/,t) 2 =t J dxf(x) 2 +M 1 (f,t) (A.15) 

with both M(f,t) and Mi(f,t) martingales. This implies that M(f,t) is Brownian motion 
in t with covariance f dxf'(x) 2 . We refer to |0| for a more complete discussion. 



APPENDIX B: 

We consider the Calogero-Sutherland model on a finite ring of length £. The density is 
p = N/t. We scale the positions with e~ l . Setting e = 2n/£, we scale back to the circle 
[0, 27t]. Then N = e~ 1 2np. Since the density scales, we set 27rp = 1 and use N instead of e. 
After these transformations the ground state ip 2 = e~^ v with the logarithmic potential 

V = -J2 log|sin(a;,-^)/2)|. (B.l) 

i<j=l 

The drift is then 

8 3 N 

%(z) = -q-PV(x) = |^cot((^ - Xi )/2) . (B.2) 

3 i=l 



For this model we show the validity of the substitution ( 5.4j) without averaging in t. 
More explicitely, we have to show that 

N N 

-^r, ; f,-)/' S . ) = — P<x*{( Xj - x l )/2)f\x 3 ) ~ ^g{ Xj ) , (B.3) 



where 



j=i 



2tt 

p 



9{x) = dycot((x- y)/2)f(y) , (B.4) 
o 

in the sense of the Cauchy principal part. Since we are on a ring, / has discrete Fourier 
coefficients fk = (l/27r) dxe lhx f(x). We assume that ^i fe i \k\ 3 fk < oo. Note that gu = 
— {f3/2)\k\fk- To show ( |B.3| ) we use a result of Johansson [^D| which states 
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and, provided h = = / , 

N N 



1 " 



Jim ((£ /(**))> = I E 1*1**/* 

j = l 8=1 ^ fc 



Let us work out the square. The first term is 

N 



which by (|B.6|) tends in the limit iV — > oo to 



§Ei*ii&i a = §£i*iUi 2 - 
^ k k 

For the second term we use -§-e~^ v = a,-e _ ^ v . Then 

JV JV AT „ 

j'=l i=l *>i=i ^ 

2 AT JV 

j=i i=i j 

By ( p.6| ) the first summand vanishes and by (P-5| ) the second summand converges 

2Y,\k\ 2 mk = -2pJ2\ k \ 3 \fk\ 2 . 

k k 
For term number three we use that 

(dP/dXidx^e-W = a^e-^ + {dcn/dxfie-W . Then 

N N 
j=l i=l 

i,j=l 1 XJ J 
Af N N 

= ^(<(£ /"(*;))(£ /''(^))> + (£r(^) 2 > + H^n^rM) 

i=i i=i i=i i 

1 fl * 

+ ]v4 <£ i sin « Xi - ^)/2)r 2 [/'(^) - z'^)] 2 ) • 
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By ( p.5| ) the first summand vanishes. The second summand converges to 

^(2 7 r)- 2 | dx J ^|sin((x-y)/2)|- 2 (f(a;)-/'(y)) 2 = ^|fc| 3 |A| 2 . (B.12) 
oo k 

Adding the three terms we conclude that the sum vanishes. 
Acknowledgements: I thank M. Prahofer for most helpful discussions. 



24 



REFERENCES 



[1] F.D.M. Haldane, General relation of correlation exponents and spectral properties of 
one-dimensional Fermi systems: application to the anisotropic S = \ Heisenberg chain, 
Phys. Rev. Lett. 45, 1358-1362 (1980) 

[2] F.D.M. Haldane, Luttinger liquid theory of one dimensional quantum fluids: I. Prop- 
erties of the Luttinger model and their extension to the general ID interacting spinless 
Fermi gas, J. Phys. C 14, 2585-2919 (1981) 

[3] S. Tomanaga, Remark on Block's method method of sound waves applied to many 
fermion problems, Progr. Theor. Phys. 5, 544-569 (1950) 

[4] J.M. Luttinger, An exactly soluble model of a many-fermion system, J. Math. Phys. 4, 
1154-1162 (1963) 

[5] K. Schonhammer, Interacting fermions in one dimension: the Tomonaga- Luttinger 
model, preprint, |cond-mat / 971033H| 



[6] D.C. Mattis, E.H. Lieb, Exact solution to a many-fermion system and its associated 
boson field, J. Math. Phys. 6, 304-312 (1965) 

[7] G. Benfatto, G. Gallavotti, A Procacci, B. Scoppola, Beta function and Schwinger 
functions for many fermion systems in one dimension. An anomaly of the Fermi surface, 
Comm. Math. Phys. 160, 93-172 (1994) 

[8] G. Benfatto, V. Mastropietro, private communication 

[9] M. den Nijs, K. Rommelse, Preroughening transitions in crystal surfaces and valence- 
bond phases in quantum spin chains, Phys. Rev. B 40, 4709-4734 (1989) 

[10] L. Balents, M. Kardar, Roughening of anisotropically reconstructed surfaces in the Hub- 
bard model, Phys. Rev. B 24, 16031-16044 (1992) 

[11] B. Sutherland, C.N. Yang, CP. Yang, Exact solution of a model of two-dimensional 

25 



ferroelectrics in an arbitray external electric field, Phys. Rev. Lett. 19, 588-591 (1967) 

[12] G. Albertini, S. Dahmen, B. Wehefritz, The free energy singularity of the asymmetric 6- 
vertex model and excitations of the asymmetric XXZ chain, Nucl. Phys B 493, 541-570 
(1997) 

[13] H. Spohn, Interacting Brownian particles: a study of Dyson's model, in: Hydro dynamic 
Behavior of Interacting Particles Systems, IMA Volumes in Mathematics and its Appli- 
cations 9, editor G. Papanicolaou. Springer Verlag, Berlin 1987 

[14] M. Prahofer, H. Spohn, An exactly solved model of three-dimensional surface growth in 
the anisotropic KPZ regime, J. Stat. Phys. 88, 999-1012 (1977) 

[15] F. Dyson, A Brownian motion model for the eigenvalues of a random matrix, J. Math. 
Phys. 3, 1191-1198 (1962) 

[16] H. Spohn, Equilibrium fluctuations for interacting Brownian particles, Comm. Math. 
Phys. 103, 1-33 (1986) 

[17] H.T. Yau, Logarithmic Sobolev inequality for lattice gases with mixing condition, Comm. 
Math. Phys. 22, 63-80 (1991) 

[18] H. Spohn, Dyson's model of interacting Brownian motions at arbitrary coupling strength, 
Markov Processes Rel. Fields 4, 649 - 662 (1998) 

[19] K. Johansson, Shape fluctuations and random matrices, preprint (1999) 

[20] K. Johansson, On the fluctuations of eigenvalues of random hermitian matrices, Duke 
Math. Journ. 91, 151-204 (1998) 

[21] H.T. Yau, S.R.S. Varadhan, Diffusive limit of lattice gas with mixing condition, preprint 
1998 

[22] A. Neumayr, W. Metzner, Reduction formula for fermion loops and density correlations 
of the ID Fermi gas, preprint (1999) 



26 



[23] C. Kipnis, C. Landim, Scaling Limits of Interacting Particle Systems. Grundlehren Vol. 
320, Springer, Berlin (1999) 



27 



